Considerable non-local electronic correlations in strongly doped Na x Co02 
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The puzzling electronic correlation effects in the sodium cobaltate system are studied by means 
of the combination of density functional theory and the rotationally invariant slave boson (RISB) 
method in a cellular-cluster approach. Realistic non-local correlations are hence described in the 
short-range regime for finite Coulomb interactions on the underlying frustrated triangular lattice. 
A local Hubbard U is sufficient to model the gross in-plane magnetic behavior with doping x, 
namely antiferromagnetic correlations at intermediate doping and the onset of ferromagnetic order 
above x>3/4 with a mixed phase for 0.62<x<3/4. Important insight is thereby provided by the 
occupations of local cluster multiplets retrieved from the RISB framework. The extended modeling 
of the x>2/3 doping regime with an additional inter-site Coulomb repulsion V on an experimentally 
verified effective kagome lattice allows to account for relevant charge-ordering physics. Therewith 
a fluctuating charge- density- wave state with small quasiparticle weight and a maximum in-plane 
magnetic susceptibility may be identified at x~3/4, just where the magnetic ordering sets in. 

PACS numbers: 71.27. +a, 71.30.+h, 71.10.Fd, 75.30.Cr 
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I. INTRODUCTION 

The investigation of frustrated spin systems plays a vi- 
tal role in condensed matter physics due to its immediate 
importance for the understanding of competing interac- 
tionsi However many realistic frustrated systems are 
insulating and may not be metalized via doping. Hence 
studying frustration effects within the metallic regime 
and possibly tuning interactions by additional doping can 
often only be realized in a pure model context. Nonethe- 
less, there is growing awareness that such effects are also 
important in correlated metals 2 -. For instance it is be- 
lieved by many researchers that for the recently discov- 
ered high-temperature superconducting iron pnictides 3 -, 
spin and/or orbital frustration is a relevant feature 4 ^ 5 -. 

Concerning the doping possibility of a strongly cor- 
related metal with manifest geometrical frustrations, 
sodium cobaltate Na x Co02 serves as an unique realis- 
tic test case. This system consists of stacked triangu- 
lar C0O2 layers, held together by Na ions inbetween. 
Nominally, the oxidation state of cobalt ranges between 
Co 4+ (3d 5 ) and Co 3+ (3<i 6 ), depending on the doping x. 
Because of the strong (t2 9 ,e g ) crystal- field (CF) splitting 
a low-spin state is believed to be realized for the Co ion, 
where x controls the final filling of the rather localized 
t2 g manifold. Due to the additional trigonal ^-internal 
aip-e^splitting, x thus formally mediates in a simplified 
view between a remaining half-filled a\ g orbital (x=0) 
and the band-insulating regime (x=l). In this respect 
however the trigonal CF splitting appears to be a crucial 
parameter. Calculations based on the local density ap- 
proximation (LDA) yield a value A(ai g — e f g )~— 0.1 eV for 
x=0.5, decreasing in magnitude away from this doping 
level^ 7 -. Yet the measured Fermi surface (FS) shows only 
a single distinct hole-like hexagonal sheet, whereas LDA 
predicts additional e' g hole pockets. Further investiga- 
tions have underlined the subtleties associated with this 
CF splitting for the cobaltate system 8 - - — . Note that a 



single hole-like FS sheet would transfer into a low-energy 
single-band modeling with a negative nearest-neighbor 
(NN) hopping t. Already on this level, an interesting 
interplay between the various degrees of freedom (spin, 
orbital and charge) is expected, where the geometrical 
frustration together with a large ratio U/W of the Hub- 
bard U and the bandwidth W should result in a complex 
behavior with x. 

Experimental investigations have indeed revealed a 
very rich phase diagram, incorporating various rather dif- 
ferent many-body phenomena and phases. A supercon- 
ducting dome (T c ~4.5 K) appears close to #=1/3 when 
intercalation with H2O is allowed for—, while a metal- 
insulator transition (MIT) is founcU^ for £=0.5 below 
^mit~51K. Above £=0.5 a disproportionation of charge 
between the Co sites takes plac o 14 i 15 and a regime of 
large thermopower is detected^ for 0.71<£<0.84. In 
view of the possible frustration effects, the varying mag- 
netic behavior is a central concern. For low doping in the 
range £<0.5 an overall Pauli-like susceptibility is founcU^ 
with however hints for two-dimensional (2D) antiferro- 
magnetic (AFM) correlation s 15 ! 17 ! 18 . Such AFM correla- 
tions are identified to increase towards x=0 for the struc- 
turally similar Li x Co02 system^. Yet isolated C0O2, 
i.e., £=0, appears to be metallic, whereby the degree of 
electronic correlations is still a matter of debat o 19 i 2Q . The 
case £=0.5 seems exceptional, with in-plane AFM order- 
in the insulating regime. For £>0.5 the system displays 
increasing magnetical response with a spin-fluctuation 
regime in the range 0.6<£<3/4, including the evolution 
to Curie- Weiss behavior 13 for 0.6<x<0.67, and the even- 
tual onset of in-plane ferromagnetic (FM) orde r 15 ! 22 — . 
The magnetic structure in the doping range 3/4<£<0.9 
with ordering temperatures Tn ^19-27K— ~ 24 i 27 is of A- 
type, i.e., the FM C0O2 layers are coupled antiferromag- 
netically. 

The fact that not the low-doping regime with the nom- 
inal Co 4+ (5=1/2) state (close to £=0) but the high- 
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doping region with the nominal Co 3+ (S=0) state (close 
to x=l) is the magnetically active part in the phase di- 
agram has motivated various theoretical efforts. Early 
LSDA computation o 31 i 32 predicted an itinerant FM solu- 
tion for all doping levels between 0.3<x<0.7, clearly in 
contradiction with experimental findings. In more recent 
calculations Johannes et al£^ were able to show that an 
A-type AFM structure is slightly favored over the FM 
state for a simple approximant to the realistic cobaltate 
crystal structure at £=2/3. Generally the different or- 
dered magnetic LSDA states are indeed rather close in 
energy 3 -^. 

Aside from the already named indications, the rel- 
evance of correlation effects beyond LDA due to the 
partially filled Co (3d) shell has been motivated by sev- 
eral experiments, e.g., from optics^, photoemission 3 - 5 - - — 
and transport^ 3 - measurements. In view of the cobaltate 
physics, Merino et discussed in detail the electron- 
doped single-band Hubbard model with NN hopping on 
the triangular lattice within single-site dynamical mean- 
field theory (DMFT). There it was concluded that for t>0 
the paramagnetic (PM) Curie- Weiss metallic phase has 
an instability towards FM order for large enough U/W. 
However such an transition is missing for the case £<0, 
where moreover the PM metallic phase displays Pauli-like 
magnetic response. In Na x Co02 the dominant NN hop- 
ping appears indeed to be negative, but relevant higher 
hopping amplitudes (at least to 3rd order) are necessary 
to bring the tight-binding (TB) modeling in good accor- 
dance with the LDA dispersio n 40 ! 41 . These higher contri- 
butions introduce a flatness to the low-energy band, giv- 
ing rise to substantial additional structure in the density 
of states (DOS) compared to the NN-TB DOS^^. This 
allows for a more interesting magnetic behavior with elec- 
tron doping in contrast to the simplified NN-TB (£<0) 
model. In fact the study of Korshunov et al^ 3 - based 
on the combination of the realistic three-band (a± g , e' g ) 
dispersion with a Hubbard on-site interaction resulted 
in (q=0)-peaks in the magnetic susceptibility at around 
Xp^0.56 (U=0 eV), x p ~0.6 (infmite-Z7 Gutzwiller) and 
Xp~0.68 (infinite- U Hubbard I). Gao et aL 44 obtained 
in similar calculations using a single-band third NN-TB 
model also within infinite- U Gutzwiller a correspond- 
ing value x p ^0.67. In a charge-selfconsistent doping- 
dependent LDA+Gutzwiller treatment with three corre- 
lated (dig, e' g ) orbitals the A- type AFM order was found 
to be stable for x>0.6 (17=3 eV) and x>0.5 (17=5 eV) by 
Wang et al. 9 , but with additional findings of FM order 
for x<0.3. Note that all these correlated studies neglect 
non-local effects of strong correlation, i.e., the self-energy 
is taken to be local and, importantly, a single correlated 
Co site in the respective planar part of the unit cell is 
assumed. Moreover the doping dependence of the sys- 
tem is simulated via simple band filling or the virtual 
crystal approximation (VCA) 4 ^. In a single-site DMFT 
study by Marianetti and Kotliar the importance of the 
Na positions for the strong correlation effects was actu- 
ally demonstrated 4 ^. 



Interestingly, another cobaltate modeling 4 ^ applying 
the t-J model with only NN-TB (t<0) hopping (and sim- 
ple doping description), solved via exact diagonalization 
on clusters as big as 18 sites (with periodic boundary 
conditions) yield Curie- Weiss behavior for x<3/4 in con- 
trast to the single-site DMFT Hubbard-model studies 3 - 9 -. 
In addition recent experimental findings 48 show evidence 
for a kagome-lattice imprint due to intricate Na ordering 
at #=2/3. These are only two hints towards the need 
for a reinforced description of non-local correlations in 
the cobaltate theory, strengthening also the original as- 
pect of the relevance of the geometric frustrations. Both 
advanced aspects, i.e., strong non-local correlations and 
a more elaborate doping dependence within a realistic 
framework were combined in Ref. 49. There the realis- 
tic LDA dispersions for a three-site Na^CoC^ in-plane 
cluster were supplemented with a finite Hubbard inter- 
action in a cellular-cluster application of the rotation- 
ally invariant slave boson (RISB) metho d 50 ! 51 . Note that 
there in the construction of the doping-dependent LDA 
Hamiltonian the respective crystal fields within the tri- 
angular cluster (modified by close-by Na ions) were not 
site-averaged as in the VCA approaches. 

In the present work we want to extent the considera- 
tions of Ref. |H by a detailed investigation of the cluster 
multiplets in the realistic three-site cluster approach to 
Na^CoC^. Galanakis et al. 52 performed a similar study 
for the doped Hubbard model on the triangular lattice 
using the non-crossing approximation (NCA) to cellular 
DMFT, but without incoporating realistic band disper- 
sions and without elucidating magnetism. 



II. THEORETICAL APPROACH 
A. Realistic modeling 

We describe the sodium cobaltate system from a 
strongly correlated perspective within a realistic frame- 
work. This however does not presume that Na^CoC^ 
may be pictured as a doped Mott insulater, a view- 
point often taken for the physics of high temperature 
superconductivity 5 - 3 -. Instead here the interplay of cor- 
relation and frustration effects in an itinerant Fermi- 
liquid (-like) state with the important proximity to the 
band-insulating (x=i) state appear to be the main char- 
acteristics. The Fermi-liquid character renders the LDA 
description as a reasonable initial starting point. Fur- 
thermore the chemical aspects due to the intriguing or- 
dering behavior of the sodium ions across the x-T phase 
diagram brings additional importance to a careful inves- 
tigation of the realistic electronic structure with doping. 
Still because of the narrow bandwidth of the low-energy 
sector, even moderate values of the local Coulomb repul- 
sion on a given Co site puts one in the regime U/W>1. 
Thus neglecting this fact may not be adequate on the 
way to a complete theory of this materials system. 

It is aimed for a non-local approach in order to al- 
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FIG. 1: (Color online) ai^-like Wannier-like functions in the 
Na x Co02 super cell for x=2/3 . Large blue (dark) balls de- 
note Na ions, light blue (light gray) balls Co ions, small red 
(gray) balls mark O ions. 



low for inter-site self-energy effects on the triangular lat- 
tice. To this, as a first approximation, the minimal 
three-site cluster with on-site Hubbard interactions is 
identified as the canonical model. For the lattice de- 
scription a cellular-cluster methodology shall serve as 
method of choice to combine the NN inter-site correla- 
tions with the overall itinerant behavior. This extended 
approach beyond the single-site description is supple- 
mented with realistic Kohn-Sham (KS) dispersions from 
LDA in a Wannier-like basis utilizing the maximally- 
localized techniqu e) 54 ! 55 . The underlying band-structure 
calculations are performed using a mixed-basis pseudopo- 
tential implementation^. This method 5 -^ utilizes norm- 
conserving pseudopotentials with a combined basis set of 
plane waves and non-overlapping localized functions. 

In the construction of the realistic model the question 
concerning the correlated subspace, i.e., the nature and 
number of included correlated orbit als, arises. In princi- 
ple the full (dig, e' g ) tig manifold appears as the suitable 
minimal choice. Yet treating a metallic state by a three- 
site cluster with altogether nine interacting orbitals and 
reduced symmetries in a large doping range is numer- 
ically very heavy for any well-suited many-body tech- 
nique. But note that the experimental FS studies reveal 
only a single a^-like band at the Fermi level and also the 
basic LDA studies predict the e^-like bands close to com- 
plete filling. In addition, in this work we are mainly in- 
terested in the magnetic behavior at larger doping where 
in any case these e^-like bands sink even deeper in en- 
ergy. Hence we reduce the correlated subspace to three 
effective a± g orbitals per Co cluster, to be determined by 
the Wannier-like construction and therefore including hy- 
bridization effects with the remaining orbitals. This ap- 
proach seems also suitable from a local-correlation point 
of view since the (nearly) filled e' g orbitals are in a singlet 
state (low-spin Co) and thus do not give rise to addi- 
tional spin scattering. In other words, double-exchange 
processes are expected to be irrelevant. At small doping 
however the inclusion of the e' g orbitals seems necessary 
because of an apparent increased relevance of the multi- 
orbital character. For the interesting magnetic properties 
at larger doping one expects the non-local correlations to 
dominate any intrasite multi-orbital correlations in this 
specific filling scenario. 

A meaningful and numerically efficient characteriza- 



tion of the complex doping behavior may be build on 
LDA calculations at selected doping levels. A super- 
cell involving the NN Co triangle of a given C0O2 layer 
serves as the base structure. By decorating this trian- 
gle with Na ions above the Co sites, i.e., usually called 
Na(l) position (see Fig. [I]), we may represent the dopings 
x={0, 1/3, 2/3, 1}. Note that in this simplified ansatz no 
sodium ions above/below oxygen positions, i.e., Na(2) 
positions, are considered. Albeit this limits the compara- 
bility to the true structure s 48 ! 58 ! 59 it allows for a discrimi- 
nation of the Co sites with x and shall be sufficient for the 
overall behavior with doping. In addition we project from 
the start the obtained complete KS Wannier Hamiltonian 
to the 2D sector and treat only the in-plane dispersion 
within a given layer. Finally the doping-dependent 3x3 
cluster band Hamiltonian is constructed through linear 
interpolation: 

H KS (K,x)= ; rHK ( f(K) + (l-a ; )H^(K) , (1) 

where x a ,Xfo are the neighboring LDA-treated dopings 
and K denotes the wave vector in the supercell for- 
malism. Note that while VCA averages in real space, 
this kind of interpolation takes place in reciprocal space. 
Thereby the local differences between the Co sites remain 
vital throughout the doping scan. Note that the 2D-only 
approach is generalized to the full 3D case within the 
effective-kagome-lattice modeling discussed in section [Vj 
The complete many-body Hamiltonian for a given x to 
be solved within the cellular-cluster approach then reads 
as follows 

H = J2 dl ia d Kj<7 + J2 H c > (2) 

Kijcr Oi 

where z,j are sites within the triangular cluster a, a is 
the spin index and is the electron annihilation (cre- 
ation) operator. The local cluster Hamiltonian H a shall 
consist of an quadratic on-site part, the local Hubbard U 
interaction term and possibly a NN inter-site Coulomb 
interaction term with amplitude V, i.e., 

H a = ^ £gjj d\ a d ja + U n^riii + j ^ n ia n ja > , 

ijcr i i^j,cra f 

(3) 

with ni a =dl a d ia . To include the Coulomb coupling be- 
tween NN clusters a mean-field (MF) decoupling of the 
inter-cluster interactions is performed, resulting in an 
embedding, which reads (dropping here the spin indices 
for convenience): 

^riirij -^(niirij) +nj(ni) - (n i )(n j )) - ^n^rij) 

(4) 

where i is a site on the examined cluster and j is a site 
on the NN cluster. In the last step leading to eq. (j4]) 
it is assumed that the average value for nj may be in- 
serted. We neglect non-density-density terms in this MF 



decoupling 6 - - since their effect should be minor. Hence 
the modified version of eq. (|3]) reads 

ija i 




The factor 1/2 in the MF term is needed to avoid the 
double counting of the inter-cluster Coulomb interaction 
when summing over a in eq. (|2]). In the actual calcu- 
lation the values (n Jcr /) are initially set and then deter- 
mined self-consistently via an outer loop to the RISB 
scheme. In a PM homogeneous case, the MF (nj a >) is in- 
dependent of j and a' and hence the contribution of that 
term corresponds to a mere chemical-potential shift. The 
justification of an (extended)-Hubbard-type Hamiltonian 
for sodium cobaltate may not straightforwardly be veri- 
fied, because additional interaction terms could arise in 
a minimal model. However we believe this approach to 
be reliable for most of the key physics in cobalt ates. Ex- 
plicit exchange terms due the fact that we integrated out 
the oxygen states are expected to be small in size and are 
very difficult to compute for metals in a rigorous man- 
ner—. 

The dispersions are provided by the KS Hamiltonian 
of eq. (pQ) (spanned on Nk points in reciprocal space) in 
the following way: 

SKij = H$ a (K, x) - ^- Y, H ™ ( K > *) > ( 6 ) 

Thus the dispersive part carries no local terms, all of 
those are transfered into H a . In the solution of the prob- 
lem one has to be aware of the breaking of translational 
invariance in the cellular cluster scheme (for reviews of 
cluster approaches to strongly correlated systems see e.g. 
Ref. j6ll463j). Since no detailed K-dependent properties 
are discussed, no periodization of the resulting cluster 
self-energy is performed in the present work. 

B. RISB method 

The full Hamiltonian from eq. (|2]) is solved employ- 
ing the rotationally invariant slave boson (RISB) for- 
malis m 5Q i 51 in the saddle-point approximation. Within 
RISB the electron operator d ia is represented as 
d ia =R[(j)]^ J ' fj a >, where R is a non-diagonal transforma- 
tion operator that relates the physical operator to the 
quasiparticle (QP) operator f ia . The transformation R 
is written in terms of the slave bosons {^An}- I n this gen- 
eralized slave-boson theory (j) carries two indices, namely 
A for the physical-electron state and n for the QP Fock 
state. It follows that the kinetic if -dependent part of eq. 



(j2j) is expressed via the QP operators with renormalized 
dispersions. The operator character of the local cluster 
part is described solely via the slave bosons^. In order 
to facilitate this operator decomposition into an QP part 
and the interacting part carried by the slave bosons, two 
constraints have to be enforced. They shall normalize 
the total boson weight to unity and ensure that the QP 
and boson contents match at every site, respectively: 

E^An = 1 , (8) 

An 

E4^(«l/Uv|n') = /l/y • (9) 

Ann' 

This selection of the physical states is imposed through 
a set of Lagrange multipliers {A}. In the mean- field ver- 
sion at saddle-point these constraints hold on average, 
with the bosons condensed to c numbers. The solved 
saddle-point equations for {</>}, {A} yield a free-energy 
F of the interacting electron system that may be used 
to investigate phase stabilities. Of course, this free en- 
ergy shall not be confused with the total free energy of 
the system that also incoporates further terms, e.g., the 
nuclei contributions. However already therefrom, many 
important insights in the possible instabilities can be re- 
trieved. Note that the K-point integration with a smear- 
ing method (Fermi, Gaussian) introduces a small effec- 
tive temperature into the system, which however may not 
easlily be interpreted as a true physical temperature due 
to the fully-condensed boson treatment in this simplest 
mean-field approach. 

The RISB physical self-energy of the d electrons at 
saddle-point is given by^ 

E d (cj) = u (1 - Z- 1 ) + [Rt^ART 1 - s (0) , (10) 

where Z=RR^ is the QP- weight matrix and A is the 
lagrange- multiplier matrix. The quantity is identi- 
cal to the one-body term in the local Hamiltonian and 
hence in the present case its elements are given by eq. (O . 
Thus Yid contains only a term linear in frequency and a 
static part. Though more approximative than involved 
methods like, e.g., quantum Monte- Carlo, which may 
handle the full frequency dependence, this is sufficient 
to describe the QP nature at low-energy. It also gives an 
idea of the effect of local excitations within a static self- 
energy. In this regard, the slave-boson amplitudes yield 
direct access to the occupation of local multiplets in the 
metallic state. The RISB mean- field approach is more- 
over rather efficient with very reliable qualitative (and 
often even good quantitative) results in most cases. Note 
that though a lattice implementation of the this method 
is used here, the technique may equally well be utilized 
as an impurity solver for sing le-site DMFTH. Due to the 
inherent local nature at saddle-point, this specific DMFT 
impurity solution is then identical to the direct lattice- 
calculation result. 

In the present case the formalism is employed within a 
cellular-cluster context, i.e., the off-diagonal components 
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of the self-energy Ed(cj) describe inter-site correlation ef- 
fects between the Co atoms in the triangular cluster. To 
obtain an intuitive interpretation for the slave-boson am- 
plitudes, we rotate both the cluster- (A) and the QP-(n) 
index of (j)An into the eigenbasis T of the isolated cluster 
problem. The transformation matrix U(T) is obtained 
via simultaneous diagonalization of the commuting set of 
operators {H ai S 2 , S z , TV}, where S is the total cluster 
spin operator, S z its z-component and TV the particle- 
number operator. Since H a is a real, symmetric matrix, 
one ends up with a representation 

(j) TT , =U ] VA (j) An U nT , (11) 

in which the slave bosons (j) rr , are directly tailored to 
the multiplet structure of the local cluster a. Note that 
the given set of quantum numbers is still not sufficient 
to classify the cluster eigenstates completely. In order 
to simplify the presentation we limit however the discus- 
sion to this definition of the T basis, which is ample for 
the physics. In the following, a specific multiplet will 
be denoted T™ n , where p describes the particle sector, 
n the energy level therein (starting with n=0 for the re- 
spective ground state) and m marks the spin state, i.e., 
singlet V, doublet 'd', triplet V and quartet 'q\ Where 
quantum numbers are stated for a given multiplet, we 
follow standard Dirac notation \E p ^ n ,S,S z ) with E the 
multiplet energy, and S,S Z total spin and z-component 
respectively. 

III. NN-TB MODEL ON THE TRIANGULAR 
LATTICE WITH LOCAL COULOMB 
INTERACTION 

To give an introduction to our approach and to touch 
base with other recent work for the Hubbard model on 
the triangular lattice, we review here results of an NN- 
TB model with t>0. Hence V is put to zero in this 
section. This model considerations shall set the stage 
for a subsequent discussion of a more realistic model of 
Na^CoC^ with higher hopping terms and, importantly, a 
NN hopping t<0. We first investigate the isolated trian- 
gular cluster to understand the structure of the multiplet 
states and then relate these findings to the physics of the 
itinerant system. 

A. Cluster states in the local limit 

A detailed study of the Hubbard model on the triangu- 
lar cluster, comparing the cases t^gO, has been presented 
in Ref. [39|. There it was shown that the tendency towards 
ferromagnetism is indeed stronger for t>0 already in the 
local quantum- chemical problem. Here we only sketch 
the t>0 case and discuss the t<0 one in line with the 
realistic sodium cobaltate model in section HVl 

Lets start with the relevant multiplet states on the 
non-interacting triangle within the different particle sec- 
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FIG. 2: (Color online) Triangular cluster spectrum of H a for 
U=0t and £>0. The height of each peak shows the respective 
multiplicity. 



K A A A 
AAA 

E 2Q ,S=0, S z =0 


Ai 1l i A 

E 40 ,S=1, S z = + 1 


'"AAA 

-A + A +A 

i + mi n 

E3„,S=1/2, S-1/2 


r d 

5,0 

4-tt t ft A 

E 5,0' S=1/2 ' S z=" 1/2 



FIG. 3: Relevant cluster multiplets for the non- interacting 
TB model. The multiplets are shown as a superposition of 
Fock states. For pictorial simplicity, numerical factors are 
suppressed. Because the factors depend on U and x, the Fock 
states are sorted after the absolute value of the factors de- 
scending from left to right and top to bottom. 



tors, obtained from a simultaneous exact diagonalization 
of {H a , S 2 , S z , N}. Note that in the present context H a 
corresponds to eq. ([3]) with e ai j=—t(l — Sij). Figure [2] 
displays the energy spectrum of the triangle in the non- 
interacting case, the decomposition of the correspond- 
ing eigenvectors into Fock states is shown in Fig. [3] for 
some relevant multiplets. For pictorial simplicity only 
the Fock states participating in a given multiplet are 
discussed, whereas the respective expansion coefficients 
are not shown, but of course are taken care of exactly 
within the calculations. In order to give an idea about 
the relative size of these coefficients, the pictorial expan- 
sion terms in Fig. [3] (and in the similar following figures) 
are ordered by the decreasing absolute value of the ex- 
pansion coefficients. 

The triangle spectrum of the NN-TB model displays 
a high degree of symmetry around zero energy in the 
non-interacting case. The lowest state in energy is a 
two-particle singlet, the highest state is given by a four- 
particle singlet. The isolated triangle with NN hopping 



6 



: 1° 


i 


i 






c 

! 4,1 






r 


5,0 

It 


(a) 

T 


-•- Op singlet 

— lp doublet 
2p singlet 

-*- 2p triplet 

— 3p doublet 






i 


i 




r 


ii 






r 


(b) 

d 

5,0 

M 


-+- 3p quartet 
-•- 4p singlet 
-*- 4p triplet 

5p doublet 
-•- 6p singlet 

■ ■ , I 





10 20 30 

E[t] 
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and U =0t is anti-symmetric under exchange of particles 
and holes, two doublets from the three-particle sector 
show up at E=±3t. The overall ground state is fol- 
lowed by the three-particle doublets, one-particle dou- 
blets and degenerated four-particle singlets and triplets 
at E=—2t. As pointed out in Ref. [39|, the three different 
spin-orientations of the triplet in the four-particle sec- 
tor are degenerated and the S z state may be changed by 
spin flips without energy cost. The resulting high sus- 
ceptibility for spin polarization in that particle sector is 
lost for the case t<0. Figure [3] shows that for t>0 the 
structure of the ground state in the four-particle sector 
is rather simple (but becomes more complex in the t<0 
case^). The simplicity of the dominant doublet Tg in 
the five-particle sector is due to the distribution of a sin- 
gle non-local hole, while the three-particle doublet 



as well as the two-particle singlet T| display a more 
complex decomposition into Fock states. 

For t>0 the electrons on the triangle can lower their en- 
ergy via hopping. The multiplet T| contains nine Fock 
states and each electron in those states has two options 
of hopping, collecting one —\t\ every time. From a chem- 
ical point of view, a bonding state can always be formed. 
In the case of the three-particle ground state T^q, the 
unpaired electron enters an anti-bonding configuration, 
thus leading to a higher state energy. 

Figure [H shows the triangular cluster spectrum in the 
interacting regime for U=7t and U=10t, the relevent 
multiplets in Fock-space decomposition are shown in 
Fig. \5\ For the investigated U values, the two-particle 
singlet remains the overall ground state, because it only 
has to pay the U penalty for double occupation three 
times out of nine. It is energetically followed by the one- 
particle doublet which does not gain that much energy 
from hopping processes because of the single-electron oc- 
cupation. But it now benefits from being free of U con- 
tributions. Hence with increasing U the overall ground 
state will eventually evolve towards this one-particle state 
via a "local Mott transition" . Note that the two-particle 
triplets are also not affected by U because they are also 
free of double occupations. The three-particle doublet 
which was following T| f° r U =0 eV, now has to pay 
the Coulomb penalty six times out of nine. On the other 
hand the three-particle doublets that show the high de- 
gree of multiplicity for the non-interacting case at zero 
energy split up for finite U. It is also interesting to re- 
mark that there is a swap in energy eigenvalues at U ~9£, 
i.e., the energy of a four-particle singlet increases faster 
in energy than the other states in this sector. Besides 
this one, no further level-crossings were observed with 
varying U. Again because of the lack of particle- hole 
symmetry, the two-particle sector is favored compared to 
the four-particle sector even at half filling. Note also that 
the four-particle triplet is lower in energy than the four- 
particle singlets. Thus in the atomic limit the triangle is 
still easy to polarize even for U>0. This holds for the 
case of three-, four- or five-particle filling, since all ener- 
getically favorable states from these particle-sectors show 
a high degree of multiplicity and have a free spin. While 
for moderate U the overall ground state T| displays a 
very complex structure it has no free spin, compared to 

1 3,0- 



B. Lattice modeling 

In the cellular framework the infinite lattice is build 
up by finite clusters, here for the triangular case given 
by the minimal triangle (see Fig. [6j). Since we are inter- 
ested in the canonical 2D lattice, the inter-cluster hop- 
ping is identical to the intra-cluster one, i.e., the disper- 
sive model is described by a NN-TB parametrization. In 
the non-interacting (Fermi gas) case the cellular formal- 
ism coincides with the standard single-site description, 
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where the bandwidth amounts to W=9t. Because of the 
absence of particle-hole symmetry the density of states 
(DOS) for smaller filling is strikingly different from the 
one at larger filling (cf. Fig. [6]). A van- Hove singular- 
ity is located at £=£ max — t (e m ax being the maximum 
single-particle energy for the single band), while such a 
singularity is missing at negative single-particle energies 
e. Note that hence importantly the chemical potential 
for the half-filled case is not placed in the middle of the 
band. 

For small Hubbard U the interacting model describes 
a Fermi liquid. We concentrate in the discussion on 
the paramagnetic regime of the half-filled as well as the 
electron-doped case. The resulting QP weight Z and 
inter-site spin correlation (Si-Sj) on the triangular cluster 
with increasing U at half filling are shown in Fig. |7{a,b). 
It is seen that the diagonal (onsite) QP weight Zi de- 
creases as U grows, until a first-order Mott transition is 
reached at U C 2~12.2t=1.36W. The off-diagonal (inter- 
site) elements Zij show opposite behaviour and increase 
towards U C 2, but are much smaller compared to Zi. The 
literature on the critical U for the Hubbard model on 
the triangular lattice, without invoking magnetic order, 
is rather large. The present critical U is in the same range 
as values from exact diagonalization of a finite cluster— 
(-12.lt) and cellular DMFT 66 (~10.5t). Figure[7]o shows 
that due to the dominant superexchange process at half 
filling the value for (Si • Sj) is always of AFM character, 
with a diverging slope at U c . 

Based on the local-cluster-limit study, the occupation 
of the extracted multiplets may now also be investigated 
in the itinerant problem. This allows for a discussion 
of the correlated metallic state for different Hubbard U 
and finite doping x from a finite-cluster real-space view- 
point. In this respect it may be observed in Fig. [8] that 
the cluster multiplet occupations are a valuable source 
of information. When discussing the condensed slave- 
boson amplitudes (j) VT , in that basis, we focus in the fol- 
lowing on the first index T and sum over the second V 
(of course, the full structure was used within the calcu- 
lations). Note that the cluster dispersion on the lattice 
leads in the present case to (weakly) nondiagonal (j) TT , 
matrices. For small U, the non-degenerate two-particle 



2.0 




FIG. 6: (Color online) 3-site cellular cluster approach to the 
Hubbard model on the isotropic triangular lattice for NN hop- 
ping t. Left: lattice tiling, right: density of states with the 
Fermi energy Ef chosen for half filling . 
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FIG. 7: (Color online) Observables for the NN-TB Hubbard 
model. U dependencies for (a) quasiparticle weight and (b) 
NN (SiSj) correlation functions at half filling. (c,d) same 
observables, but behavior with doping. 
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FIG. 8: (Color online) Relevant slave-boson amplitudes (j> 
squared with increasing U for the cluster NN-TB model. The 
amplitudes are labeled with a single r™ n index (see text). The 
notation {z, . . . , j} implies that the amplitudes for all states 
from i to j are of equal weight. 



singlet T^Q has the highest occupation at half filling. 
However note that because of degeneracies the occupa- 
tion of several other multiplets is identical, also allowing 
in principle for increased fluctuations among them. Thus 
the total occupation within such a class of multiplets, i.e., 
{T\ o, . . . , T\ 2 } may be even higher. When normalized to 
single-multiplet occupations, the three-particle doublet 
states |^3 o, V 2 ? ±V 2 ) dominate for U>9t. Close to the 
Mott transition the paramagnetic system can be found 
nearly exclusively fluctuating between those two doublet 
states. 

Figure [9] displays the slave-boson amplitudes with elec- 
tron doping x up to an investigated regime of £=0.425 in 
the non-interacting case and at fixed U=10t<U c . There 
the U=0t case shows a simple crossover from the dom- 
inant r| o a t half filling to the six-particle singlet Tg 
at larger filling (with statistical occupation of the re- 
maining particle sectors). On the other hand for large 
U the four-particle triplet T\ | 2 y becomes strongest 

above x^0.2. Note that the maximum of | 2 || 2 is 
reached around _c=l/3 which is in the very same range as 
the onset of ferromagnetism in single-site DMFT for the 
given U value 39 . The simple Stoner criterium is known 



8 



0.15- 



0.05 



0.00 



non-inter act. 

■" (|) d 5,{0,..,3} 





interacting: 




— 4»\o 








— $\{0 t ..,2} 




^4, {3,4} 




— <l) d 5,{0,..,3} 




— ^6,0 




FIG. 9: (Color online) Doping-dependent slave-boson ampli- 
tudes in the NN-TB model for U=0t and U=10t. 
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TABLE I: Average hopping integrals for the Na x Co02 clus- 
ter modeling from the effective a\ g Wannier-like construction 
shown here up to fourth-nearest neighbor (cf. Fig[BJ). The on- 
site level difference A^Na denotes the Na-induced level lower- 
ing for Co sites with sodium on top. Note that in the actual 
calculations the full 2D hopping matrices, i.e., not averaged 
and not restricted to 4th NN, were employed. The inter-layer 
hopping t± is shown for comparison and did not enter the 
RISB calculations discussed in this section. Values are given 
in meV. 



to fail in this context^, e.g. for x~l/3 the critical £7™ 
for ferromagnetism on the triangular lattice should be 
around At. Thus the cluster viewpoint provided by the 
multiplet occupations in the RISB formalism, incorpo- 
rating both, band dispersions and exchange mechanisms 
beyond Hartree-Fock, may serve as a good estimate for 
the magnetic ordering tendencies on the lattice. In this 
respect it is obvious that the large weight of the four- 
particle triplet renders the system rather susceptible to 
ferromagnetism by aligning the net spin of each basic 
triangle on the lattice. The dominance of the Y\ ^ 2 ^ 
state around x~l/3 is also in accordance with the re- 
cent findings of Galanakis et al£l in their cellular-DMFT 
study using an NCA impurity solver. 



IV. STRONGLY DOPED SODIUM COBALTATE 
WITH LOCAL COULOMB INTERACTION 

This section deals with a realistic study of Na cc Co02 
over a rather wide doping range based on projected LDA 
dispersions and explicit local on-site Coulomb interac- 
tions. Thus as in the previous section, V=0 eV holds 
here. We will show that such local Hubbard interactions 
within a cluster (or multi-site) modeling are sufficient to 
deliver a good description of the in-plane magnetic be- 
havior with doping. Coming from lower x, the PM phase 
with AFM character as well as the onset of FM order 
at x~3/4 results from this approach. A PM-FM mixed 
phase inbetween these regimes in the area 0.62<x<3/4 is 
obtained, signaling a first-order scenario for this transi- 
tion and the relevance of disorder /segregation effects for 
a deeper understanding of sodium cobalt ate at large x. 

As described in section Hi Al the realistic modeling of 
Na^CoC^ is based on the construction of an effective a± g 
minimal cluster model. It is derived from the LDA de- 
scription of chosen reference configurations (serving as 
simple approximants to the more complex true struc- 
tures) for the dopings #={0,1/3,2/3,1}. The Wannier- 
like Hamiltonians are projected onto an effective single 



layer, i.e., the RISB method is faced with a 2D dis- 
persion. The averaged in-plane hopping amplitudes ob- 
tained from the respective Wannier-like construction are 
given in Tab. HI Compared to the pure model case in the 
last sections, the NN t is now negative. It is moreover ev- 
ident that the absolute value of t decreases with increas- 
ing doping x, in the end loosing about half its magnitude 
towards x=l. This fact already points to a diminishing 
superexchange process at larger x. The higher hopping 
integrals (shown here up to fourth-NN) decay quickly 
with distance and also decrease in magnitude with dop- 
ing. Note that the on-site levels are lowered by putting 
sodium on top. This can easily be understood by the at- 
tractive Coulomb potential introduced by the Na + ions. 
In Tab. U the NN inter-layer hoppings t± are also pro- 
vided. They show non-montonic behavior with x and an 
expected slight increase of the three-dimensional (3D) 
character through the ratio t±/t with doping. However 
these and higher inter-layer hoppings are not taken into 
account in the lattice studies presented in section IIVBI 

A. Cluster states in the local limit 

Similar to the strict model case, we begin with the 
study of the states in the single-cluster limit. Here the 
non-interacting part of the Hamiltonian is given by the 
(0, 0) column of Tab. [J (note that of course just the NN t 
enters in that limit). Thus not only t<0 now holds, but 
also the symmetry of the cluster Hamiltonian is lowered 
at finite doping because of the inequivalent Co sites. Our 
aim is not to discuss the energy spectrum in the cluster 
limit at each doping level in full detail. Only the ma- 
jor differences to the high-symmetry (t>0) model case 
treated before shall be elucidated. 

The spectrum of the cluster eigenstates at x=2/3 for 
U =0 eV is shown in Fig. [101 The loss of degeneracy due 
to the lower symmetry is obvious. In general, because of 
the sign change in t, the spectrum is more or less mir- 
rored compared to the (t>0) model case. While for the 
latter, low-particle sectors where energetically favoured, 
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here the six-particle singlet at ~-3.37 eV is the lowest 
state in energy. Thereafter five-particle doublets follow 
in energy, which will turn out to be important in the itin- 
erant lattice problem discussed in the next section. As an 
additional generic difference to t>0 the particle sectors 
tend to group, larger fillings at lower energies due to the 
reduced number of possible hopping processes which are 
now penalized by energy cost. On the contrary, for the 
(£>0) case the particle sectors are distributed over the 
whole range of the spectrum for U=0 eV (cf. Fig. [2]). 

Focusing on the four-particle sector which turned out 
to be important in the doped (t>0) case, it is observable 
that now a singlet is forming the lowest energy stated. 
Remember that the onset of FM order for t>0 from the 
local viewpoint may be associated with the dominance 
of the ground-state triplet. Moreover the Fock-space 
decomposition of the four-particle ground state, shown 
in Fig. [TTJ is far more complex for the singlet than for 
the former triplet (see Fig. [3]). The ground-state singlet 
cannot be polarized. To get a polarization, the cluster 
with four electrons would have to be excited, e.g., to the 
named triplet state. But even in the non-interacting case 
this would cost an energy of about 0.3eV. 

The emerging cluster spectrum in the interacting case 
(U=b eV) for the Hamiltonian associated with the fill- 
ing x=2/3 is shown in Fig. [12j Now the lowest energy 
state stems from the three-particle sector and its highly 
degenerated doublets and quadruplets. It is followed by 
a number of two-particle states and the overall tendency 
to group states is enforced. Of course by including a 
Hubbard J7=5eV the many-particle sectors become en- 
ergetically also costly compared to the non-interacting 
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FIG. 12: (Color online) Same as Fig. [TO] but for U=5 eV. 
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FIG. 13: Relevant cluster multiplets for the ferromagnetic 
LDA-derived model. 
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FIG. 10: (Color online) Energy spectrum of the cluster limit 
of the LDA-derived model Na x Co0 2 for x=2/3 and U=0 eV. 
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FIG. 11: Relevant cluster multiplets for the paramagnetic 
LDA-derived model. 



case. In the four-particle sector, the singlet is still low- 
est in energy, but the triplet draws near, the energy gap 
becomes smaller by about a factor of three. That means 
the interacting Na^CoC^ triangle with a filling of four 
electrons is still less polarizable in the local cluster limit 
than the triangle with £>0, but the local Coulomb repul- 
sion naturally tends to support tendencies towards spin 
polarization. As shall be seen in the itinerant problem, 
the following three particle doublets and other singlets 
can be neglected. Next in energy and easy to polarize 
due to a single free spin on the cluster is the degenerated 
five-particle doublet. As shown in Fig. [13] it resembles 
the structure of the easy to polarize four-particle triplet 
from the (t>0) case. In other words, while from this 
local viewpoint the (t>0) case exhibits susceptibility to- 
wards FM order at smaller doping associated with the 
four-particle triplet, the Na^CoC^ (t<0) case shows such 
tendencies only at higher doping connected to the five- 
particle doublet. 
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FIG. 14: (Color online) QP spectral function of PM Na*Co02 
at x=l/3 and x=2/3. The splitting into three parts, respec- 
tively, is due to the cluster description (see text). 
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FIG. 15: (Color online) Doping-dependent quasiparticle 
weight (a) and local self-energy (b) for PM NaccCoCb, given 
in the eigenbasis of the QP operator i.e., the orbital 

density matrix. Close to x=2/3 the red/grey orbital, cor- 
responding to the strong filling-dependent state in the inset 
O0ci , also red/grey) shows a stronger renormalisation than 
the other two (blue/dark). 



B. Lattice modeling 

The central aim of this subsection is to investigate the 
in-plane magnetic behavior of sodium cobaltate close to 
zero temperature with doping x. It is known that even- 
tually the C0O2 planes couple antiferromagnetically (A- 
type AFM order). Here we want to address the question 
why Na^CoC^ displays different in-plane magnetic cor- 
relations with doping and finally FM magnetic order for 
3/4<x<0.9. In this respect the mean- field RISB method 
allows for the possibility to stabilize paramagnetic as well 
as symmetry-broken ferromagnetic phases for the realis- 
tic Hamiltonian (eq. ([2])) 49 . For instance, a FM phase 
may be set up via initial spin-different lagrange multi- 
pliers and this lowered symmetry is self-consistently re- 
produced at saddle-point convergence if such a phase is 
(locally) stable. U =5 eV was set for the interacting case 
and we concentrate on dopings #>l/3 where the a± g - 
like modeling should be well justified. Note that now of 
course also the LDA-derived hoppings beyond NN dis- 
tance come into play, thus a one-to-one comparison with 
the NN-only model has to be taken with caution. 

The QP spectral function in Fig. [HI shows the decom- 
position into three parts due to the cluster description. 
The different parts represent the K-integrated spectral 



weight of the resulting super cell bands. They naturally 
add up to the known ai^-like cobaltate dispersion^^ in 
the non- interacting case. Note that we only consider the 
2D dispersions in K-space. As expected from the doping 
levels, the overall renormalization is weakened at £=2/3 
compared to the £=1/3 case. In the higher doping regime 
the spectral weight strongly increases close to the addi- 
tional van-Hove singularity, which could point to the rele- 
vance of Stoner-like mechanisms for the onset of in-plane 
FM order. 

The variation of the QP weight Z and the static self- 
energy part E s (see eq. ([TO]) ) with x is plotted in Fig. [T5l 
We display the quantities in the eigenbasis of the cluster 
orbital density matrix (ft f) of the correlated problem in 
order to make the point that the usual site-basis may not 
always be the most adequate one. Note that resulting 
eigenstates do not directly correspond to the spectral- 
weight parts from Fig. [HJ since the latter stem from a 
diagonalization at each K-point. It is seen that two of the 
eigenstates of (P f) , ^2,3? are highly occupied whereas 
the filling of the remaining effective orbital ip Cl shows 
strong doping dependence. Z(x) and T ls {x) exhibit the 
strongest differences among those states around x~2/3, 
pointing towards dominant inter-site correlations. Rep- 
resented in the site basis, the splitting of the diagonal el- 
ements of Z(x) and T, s (x) is diminished. However on the 
present modeling level the discrepancy between the effec- 
tive orbit als, aside from the filling behavior, is not signifi- 
cant. The overall degree of correlation with Z~0.82-0.86 
in this designated doping regime is still modest. We will 
elaborate on this point in section [Vj 

Figure [16] shows the pair-averaged spin-correlation 
functions (SiSj) between Co sites ij on the 3-site cobal- 
tate cluster. Already the non-interacting treatment 
shows substantial AFM correlations for small doping in 
the PM phase, which are strongly reinforced by the Hub- 
bard U. The latter is easily understandable from the rele- 
vant superexchange mechanism close to half filling (x=0). 
But on the other hand interestingly the AFM correlations 
are more significantly suppressed for finite U at large 
x. At x~2/3 there is a clear crossover in these short- 
range correlations visible between interacting and non- 
interacting case. Though a generally diminished {SiSj) 
amplitude at larger doping is reasonable from the sole 
reduction of the local magnetic moments, this crossover 
remains remarkable. Hence the spin correlations in the 
correlated regime show highly non-trivial behavior with 
doping. In addition, at x^O.61 indeed a long-range or- 
dered FM metal may be stabilized. The metallic cluster 
magnetization of that phase exhibits linearly decreasing 
behavior with doping, with subtle convergence proper- 
ties for x>0.9 concerning the K-point integration with 
smearing methods. 

Figures [T71 and [TBI show the relevant slave-boson ampli- 
tudes for the competing phases with doping x. In Fig.fTTl 
only the paramagnetic phase is investigated. Since the 
system is in a high-doping regime, multiplets from the 
four-, five- and six-particle sector play the dominant 
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FIG. 16: (Color online) Doping-dependent spin correla- FIG. 17: (Color online) Doping-dependent relevant slave- 
tion function on the triangular cluster for the PM interact- boson amplitudes for the interacting (£7=5 eV) and non- 
ing (blue/dark full) and non-interacting (blue/dark dashed) interacting PM phase within the LDA-derived Na x Co0 2 
phase as well as the FM phase (red/grey full). Inset: averaged model. Inset: Comparing interacting cases with complete dis- 
cluster magnetization with x. persion (full lines) and only NN hopping (dots) around x-0.7. 



role. Starting below x~0.4, the four-particle singlet is 
strongest, in line with the ongoing NN-AFM correlations 
there. At x~0.52 the five-particle doublets however over- 
run the four-particle-singlet weight and reach a maxi- 
mum at around £=0.69. Hence at x~2/3, the probabil- 
ity to find the system in one of the two doublet states 
is about 25% each. As the doping increases even fur- 
ther the six-particle singlet, corresponding to completely 
filled cluster becomes of course dominant. As expected, 
the six-particle singlet is much stronger occupied in the 
non-interacting case for constant x. However also the 
occupation of the doublet is smaller for U=0 eV, while 
the four-particle singlet gains weight in the doping range 
£>0.44. This corresponds to an enhanced establishment 
of a net spin on the triangular cluster with sizeable Hub- 
bard U, which gives the trend towards ferromagnetic be- 
havior on the lattice. At small doping below x^O.44, the 
singlet is favored by interaction, underlining the AFM 
tendencies. It is interesting to note that the dominant 
maximum of the five-particle doublet around x^0.7 only 
little depends on the influence of the hopping terms be- 
yond NN. Albeit the latter are responsible for the van- 
Hove singularity close to the upper band edge for £<0, a 
RISB calculation with only NN t (fixed here to the value 
£=—135 meV) results mainly in the enhancement of the 
four-particle singlet weight (see inset of Fig. [T7|) . Thus 
the more distant hoppings indeed strengthen the trend 
towards in-plane ferromagnetism, but it is questionable 
if this is the only explanation for the eventual FM order. 

The occupations of the slave-boson states in the PM 
and the FM phase are compared in Fig. [181 In the fer- 
romagnetic phase, the degeneracy of the dominant five- 
particle doublets is lifted with the favored S z =l/2 one 
showing a prominent maximum at x^0.68. The proba- 
bility to find the system in that state amounts to roughly 
50%. Note that the four-particle ground-state singlet 
is not occupied at all in the FM phase, instead triplet 
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FIG. 18: (Color online) As Fig. El but now comparing inter- 
acting PM and FM phases. 



states with S z =+1 from that particle sector display some 
weight. 

In order to map out the actual stability ranges of the 
PM and FM phase, the doping-dependent free energy 
has to be studied. Accordingly Fig. [T9l shows the respec- 
tive curves obtained from the saddle-point solution of the 
RISB free-energy functional 84 . Be aware that our present 
discussion is limited to the 2D limit of Na^CoC^ and that 
very general stability statements can of course only be 
made for the full 3D case 4 ^. Note also that a full FM 
order of the 3D lattice can truly not be stabilized in the 
RISB treament of our model cobalt ate. One may observe 
in Fig.[T9lthat the FM phase is only (meta)stable around 
x~2/3. Above x^0.7 it becomes locally stable, but as 
revealed by the tie-line construction, globally a mixed 
PM-FM phase is the thermodynamic stable solution in 
the range 0.62<x<3/4. This first-order scenario for the 
onset of magnetic order is substantiated by the experi- 
mental findings of A-type AFM ordered phases starting 
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FIG. 19: (Color online) RISB free energy for the PM 
(blue/dark) and FM (red/grey) phase. There is a region of 
coexistence ranging from x a ~0.62 to x^3/4 obtained from 
the tie-line construction. 



at x~3/4 with a small but still sizeable finite magnetic 
moment (in line with our theoretical values). Moreover 
the peculiar behavior in the given doping range, espe- 
cially around x~0.7, was elucidated in several experi- 
mental works^— . The RISB FM free-energy curve dis- 
plays a minimum close to x>0.9. While in experiment 
the magnetic order is lost for x>0.9, our mean- field so- 
lution in the 2D limit still allows for in-plane FM order 
up to the band-insulating regime. But the obtained nu- 
merical solutions in that doping regime are rather fragile 
and temperature effects (somehow mimiced in the present 
RISB implementation via the smearing in the K-point in- 
tegration) are very effective in easily destroying a simple 
thermodynamic phase behavior. Remember that exper- 
imentally close to x=l strong tendencies towards phase 
separation are seen (e.g. Ref. lLU and references therein), 
which may also be detected theoretically 4 ^. 

Since the employed cellular-cluster calculations are nu- 
merically rather costly, the question arises if the ther- 
modynamic behavior from Fig. [19] may also be retrieved 
from a single-site point of view. Therefore we studied 
two alternative approaches, namely the multi-single-site 
solution of our established cobaltate model and a true 
single-site model derived from the LDA band structure 
of the single Co site unit cell of NaCo02- Note that 
the former model corresponds to the neglect of inter- 
site self-energies, however keeping site-dependent self- 
energies in the cluster unit cell. The calculations revealed 
that the multi-single-site approach yields nearly identical 
free-energy curves in the relevant doping regime, with de- 
viations becoming sizeable only at smaller x closer to the 
half-filled regime. Still, quantitatively the cluster free- 
energy remains always lower by a few ten meV. On the 
other hand the true single-site model showed similar ten- 



dencies, i.e., PM phase at intermediate and FM phase at 
high doping, but resulted in rather different phase bound- 
aries not matchable with experimental findings. Hence 
inter-site self-energies seem not essential for the overall 
integrated thermodynamic behavior on the present level 
of modeling (but surely are for the detailed electronic 
signatures). Yet a non-local description in the sense of 
site-dependent self-energies in enlarged supercell descrip- 
tions is believed to be crucial for the energetics of sodium 
cobaltate. 



V. STRONGLY DOPED SODIUM COBALTATE 
WITH EXTENDEND COULOMB 
INTERACTIONS 

So far the modeling was concentrated on the generic 
doping dependence of the correlation effects within a 
single-layer approximation to the range 1/3<x<1. In 
this last section we want to elaborate on the modeling 
for the doping regime around x~2/3, which appears dis- 
tinguished because of its strong thermoelectric repsonse 
as well as its magnetic behavior in line with the appar- 
ent phase competitions discussed in section IIVI Hence 
a more refined structural modeling taking also into ac- 
count the observed charge ordering phenomena for x> 0.5 
should shed more light upon the physics in this region 
of the phase diagram. Recent experimental findings re- 
vealed the possibility of a kagome-lattice imprint in the 
C0O2 layers at x=2/3 due to Na-ordering driven charge 
order— and we want to build up a more subtle realistic 
cobaltate model on this proposal. To include the effect 
of a charge-ordered background the inter-site Coulomb 
interaction V is now vital. As will be shown, therewith 
one is also in the position to finally describe the strong 
mass-renormalization effects in this high-doping regime, 
which were out of reach in the local- U only modeling 
(cf. Fig. H5]). 

Figure [20] shows the Na-decorated C0O2 layer of an ap- 
proximant to the rather c omp lex true Na2/ 3 Co02 lattice 
structure detected in Ref. |48| (incorporating 88 atoms in 
the unit cell). The former structure has four in-plane 
Co sites in the unit cell, where one of the sites (Co(l)) 
has a Na ion on top, i.e., in the Na(l) position, and the 
other sodium ions occupy Na(2) positions close to the 
remaining Co(2) sites. Thereby the approximant is a 
simple approach to mimic the key property of the orig- 
inal structure, namely a kagome-lattice imprint if one 
assumes the Co(l) sites to be formally in the charge- 
blocked Co 3+ state 48 . Albeit a simple LDA calculation 
surely does not yield a strong charge-ordered pattern, 
we proceed in constructing an effective kagome lattice 
model for Na 2 / 3 Co02 by projecting the resulting band 
structure onto ai^-like Wannier functions centered only 
at the Co(2) sites. Since in addition we want to incor- 
porate now also the full staggered 3D character we use 
the Co(2) sites in two adjacent C0O2 layers. Hence we 
end up with a six-dimensional model where, respectively, 
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FIG. 20: (Color online) Top view on the Na decorated C0O2 
layer of the utilized crystal approximant to the experimentally 
detected kagome- lattice imprint^. Sodium atoms are marked 
blue/dark, cobalt atoms grey and oxygen atoms red/light grey. 
The resulting effective kagome lattice is sketched via dark 
lines. 



three Wannier functions are again located on an in-plane 
triangle cluster as a basis for an effective kagome lattice 
(cf. Fig. [20]) . Note that this is a simple approximation to 
an actual blocking of the Co(l) sites through projecting 
out the explicit hoppings degrees of freedom associated 
with orbitals located on those sites. From this derivation 
of the realistic Hamiltonian the NN hopping amounts to 
t=-152 meV, i.e. is slightly larger than for the initial 
triangular lattice model at £=2/3 (cf. Tab.!]). The NN 
inter-layer hopping amounts now to t±=9 meV, i.e., is 
somewhat smaller than in the pure triangular case. Thus 
the 2D character, again loosely defined via t±/t, is rein- 
forced for the effective kagome lattice. 

In the following the emphasis is still on the physics 
within a single C0O2 layer and in the PM phase. How- 
ever the description within the full 3D structure, allow- 
ing here for inter-layer resolution, is capable of reveal- 
ing not only the experimentally-observed bi-layer split- 
ting in the correlated regime. In addition we now easily 
may stabilize the A-type AFM order in the RISB calcu- 
lations. While V is varied to look for the influence of the 
inter-site Coulomb term, the Hubbard U remains fixed to 
U =5 eV. Since the physics is now studied on an effective 
lattice where the supposingly blocked Co 3+ states have 
been projected out, the nominal filling on that lattice is 
of course not identical to the original doping level x. Its 
easily shown that the letter value is related to the effec- 
tive doping x e ff for the given kagome problem through 
x ei ?=(4x - l)/3. 

The liberty is taken to allow for a stability of the effec- 
tive kagome structure also for dopings somewhat larger 
than £=2/3. Since it is known that the charge ordering 
persists in that region, the use of the (x=2/3)-structure 
with adjusted doping level may serve as a suitable first 
approximation. In Fig. [21] the orbital filling and the 
QP weight Z of the eigenstates of (f^f) for the dopings 
£=2/3,3/4,4/5, i.e. x e ff=5/9, 2/3, 11/15, are displayed. 
The transformation to the given eigenbasis becomes now 
vital since its is easily observed that the inter-site V has a 
rather strong influence on the observables. For x e s=2/3 
the effective kagome lattice is in the ordering regime of a 
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FIG. 21: (Color online) Occupations (a) and QP weight (b) 
for the extendend cobaltate Hubbard model on the effective 
kagome lattice with increasing V in the eigenbasis of the clus- 
ter orbital density matrix. 
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FIG. 22: (Color online) Inter-site spin correlations on the 
triangular cluster of the the effective kagome lattice for the 
three different doping levels. 



charge-density- wave (CDW) instability on a lattice with 
triangular coordinatio n 60 ^ 2 " — . In the present study, this 
instability is not driven by a charge localization on the 
lattice sites i, but in the eigenstates c of the triangular 
cluster. The strongly renormalized state c\ in Fig. |2T1 
corresponds locally to a symmetrical cluster eigenstate 
of the form i/j c =A Y^=i (note that t<0 inverts the 
bonding-antibonding energy hierachy), where the ampli- 
tude A is identical for all a^-like site-centered states and 
independent of V. The remaining two cluster eigenstates 
are degenerate with non-trivial linear combination of the 
site-centered states. Note that this result may serve as a 
natural explanation for the findings of small quasiparticle 
weights in the high doping regime^— . 

Thus for large V and x e ff~2/3 the system is unsta- 
ble against a resonating- valence-bond- (RVB) like CDW 
phase. Because of the cluster choice within the cellu- 
lar description, that phase still breaks translational in- 
variance. At the actual x=2/3 level, the discrimination 
between the quasiparticle weight of the respective eigen- 
states is still visible, yet the strong signature is lost. 
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FIG. 23: (Color online) Slave-boson amplitudes of selected 
multiplet states on the triangular cluster of the effective 
kagome lattice for the three different doping levels. 



The critical value V c ~0A eV for the metal-insulator 
transition at commensurate filling is perfectly reasonable 
for the cobaltate system. Concerning the in-plane spin 
correlations, Fig. [22] remarkably shows that the inter- 
site V has a qualitatively different effect on the vari- 
ous dopings. While V for x eff =5/9 (x=2/3) strength- 
ens the AFM tendencies, at x e ^=2/3 (x=3/4) the AFM 
character is substantially weakened, with a clear trend 
towards FM correlations. For # e ff =11/15 (x=4/5) the 
NN-AFM correlations are already rather weak in the PM 
phase. Interestingly, x e ff~2/3 just corresponds to the 
actual doping value x=3/4 of the full triangular lattice 
where the onset of in-plane FM order is found experimen- 
tally. In this context it is essential to note that there are 
very strong hints that the maximum experimental Curie- 
Weiss susceptibility is indeed not at £=2/3 but closer to 
the x=3/4 doping leve l 13 i 78 . Hence besides the general 
thermodynamic considerations presented in the former 
section, the explicit inclusion of realistic charge ordering 
effects also points to the onset of in-plane FM order at 
x~3/4 with a unique underlying electronic phase. 

Considering once again the slave-boson amplitudes of 
the local cluster multiplets for the investigated doping 
levels, Fig. [231 renders it obvious that for x=3/4 the five- 
particle doublets are strongly favored by V. On the other 
hand for x=2/3 that multiplet, though still dominating, 
has to cope with a still sizeable occupation of the four- 
particle singlet state. Therefrom, the enhanced AFM 
tendencies in the latter instance compared to the V=0 eV 
case are understandable from the local cluster perspec- 
tive. This four-particle singlet state becomes even more 
suppressed for x=4/5, where experimentally already the 
in-plane FM phase is established. 



VI. CONCLUSIONS 

In this paper we employed the RISB mean-field formal- 
ism to approach the sodium cobaltate system at larger 



doping by combining realistic LDA dispersions with 
tailored many-body model Hamiltonians in a cellular- 
cluster scheme. Due to the intriguing sodium order with 
doping, a realistic, numerically feasible modeling for a 
large part of the phase diagram asks for a simple, but 
still non-trivial approximate treatment of the structural 
details. We showed that non-local correlation effects (at 
least in a short-range regime) are indispensable to ac- 
count for the phase behavior in the magnetically active 
region x>0.5. Site-dependent self-energies, triggered by 
the sodium ordering, are at minimum needed to allow for 
reasonable insight into the physics of Na x Co02. 

Emphasis was taken to study the interplay between 
band-structure effects and local cluster correlations in 
the investigation of the highlighting magnetic properties. 
It became clear that though the van-Hove singularity at 
the upper band edge (due to long-range hopping terms) 
is a vital ingredient, the apparent many-body correla- 
tions on the frustrated lattice are an additional relevant 
source for the enhanced susceptibilities. In this respect, 
the slave-boson amplitudes for the triangular cluster mul- 
tiplets proved to be valuable observables to substantiate 
the latter statement. For instance, FM tendencies for 
Na ;E Co02 (t<0) close to x=2/3 in contrast to the strong 
FM susceptibility for t>0 close to x=l/3 may directly 
be connected to the local triangular cluster correlations. 
The relevance of five-particle doublet cluster states in the 
neighborhood of x=2/3 should motivate a different view- 
point on the basic entities of the cobaltates. Namely, ap- 
proaching the magnetic properties in this doping regime 
from a larger (not necessarily triangular) cluster perspec- 
tive appears superior to the incommensurate single-site 
starting point. 

Aside from the simplified Stoner approach there are 
various more specialized considerations concerning ferro- 
magnetism in Hubbard(-like) models, e.g., ideas based on 
Nagaoka's theorem 7 ^, Kanamori's criterion^ in the low- 
density regime or ring-exchange mechanisms relevant for 
triangular plaquette o 81 i 82 . A very direct connection to 
such model ideas is complicated because of the difficult 
structural facts, however it may well be that certain as- 
pects therefrom also apply to sodium cobaltates^. 

Our calculations revealed that the doping region 
0.6<x<3/4is highly susceptible for magnetic phase sepa- 
ration, pointing towards a first-order-transition scenario 
for the onset of the A-type AFM order. These theoret- 
ical findings are in accordance with many experimental 
studies in that region 67 70 , which, e.g., revealed hints for 
disorder-induced non-Fermi-liquid behavior with possi- 
bly underlying quantum Griffiths phases 7 ^. The region 
close to the band-insulating x=l limit appears interest- 
ing for several reasons, since it is expected that a minimal 
hole doping of that stated will result in nontrivial charge 
susceptibilities compared to standard weakly correlated 
systems. However the simple structural model prohibits 
at this level detailed deeper investigations. 

The final study on the ambitious modeling of dopings 
close to x=2/3 based on the more elaborate effective 
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kagome lattice within a charge-ordered background seem 
rather fruitful concerning the understanding of the inter- 
play between the charge and spin degrees of freedom. It 
showed that the charge ordering may be efficient in sin- 
gling out distinct phases that can give rise to enhanced 
FM tendencies at specific dopings. This viewpoint deliv- 
ers a rather natural possible explanation not only for the 
onset of FM order above #=2/3, but provides a clear-cut 
interpretation for the strong correlations in that doping 
regime. The blocking of some of the charge not solely 
reduces the effective hopping. It also gives rise to size- 
able inter-site Coulomb interactions that strongly renor- 
malize the remaining half-filled effective state emerging 
from the neigborhood to an fluctuating CDW instabil- 
ity on the still accessible lattice sites. This general idea 
carries over to the specific kagome pattern and may be 
generalized to the even higher doping regions by think- 
ing of larger clusters with more blocked sites. Future 
work shall also be focused on the understanding of the 
spectral properties of such phases, probably also relevant 
in the misfit cobalt ates, in order to touch base with al- 
ready available angle-resolved photoemission (ARPES) 
datsb^^r— . In this respect, magnetic and ARPES 
measurements already suggested a fluctuating character 



of a hidden collective CDW-like instabilit y 23 ! 75 , possibly 
driven by long-range Coulomb interactions. 

In summary, because of its compelling physical rich- 
ness the sodium cobaltate system serves as a demanding 
testing ground for the nowadays available realistic many- 
body tools. Non-local correlations effects naturally show 
up in Na^Co02 and this shall stimulate further work on 
schemes to handle those as well as on many-body Hamil- 
tonians beyond the standard Hubbard model. 
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